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Abstract 


The  decay  of  the  solar  magnetic  field  nay  be  computed  frcm 
a  variational  principle  based  on  Lena's  lav;:  i.e.,  since  the  fields 
resist  any  ctianges,  they  adjvtst  themsel-'.-es  to  make  the  decay  time  a 
maxijniun,  A  variational  calculation  of  the  decay  time  associated  with 
one  mode  of  the  solar  magnetic  field  has  been  made  and  compared  v.dth 
an  exact  value  obtained  by  Wrubel.  To  the  accuracy  of  the  calculation, 
a  simple,  practically  constant,  trial  function,  plus  one  iteration,  re- 
produces the  exact  result  identically. 
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1,  Introduction 

The  decay  time  associated  with  the  solar  magnetic  field  was  investi- 
gated some  years  ago  by  Cowling'-  ■',  who  used  a  variational  method  to  calculate 

[21 
the  decay  time  associated  with  the  lowest  mode.  More  recently,  Wrubel'--'  com- 
puted decay  times  for  several  higher  modes.  Although  Wrubel  analyzed  the 
higher  modes  by  integrating  the  differential  equation  of  the  problem  numeri- 
callyj  Cowling's  use  of  the  variational  method  may  be  ectended  to  apply  to  the 
liigher  modes  also,  and  it  is  this  aspect  of  the  problem  that  will  concern  us 
here, 

A  general  discussion  of  the  ^jse   of  variational  techniques  for  the 
calculation  of  decay  times  xras  given  in  a  previous  report  ^•-',  This  included 
both  a  set  of  variational  expressions  for  computing  decay  times,  and  a  vector 
iterational  procedxire  for  systematically  improving  the  trial  functions.  It 
was  also  pointed  out  that  the  stationary  property  of  the  variational  principles 
was  the  mathematical  expression  of  the  physical  fact  that,  since  the  actual 
fields  resist  any  change,  they  adjust  themselves  to  make  the  decay  time  a 
maximum  ( Lenz '  s  lav?) . 

In  the  report  mentioned'-^-',  the  aim  was  to  present  an  exhaustive 
treatment  of  the  various  variational  principles  encountered  in  decay  problems; 
the  discussion  was  therefore  kept  general.  Here,  in  contrast,  we  vn.sh  to  dis- 
cuss  only  the  specific  problem  considered  by  Wrubel*-  ■• ,  A  single  variational 
principle  is  used  and  detailed  numerical  results  are  obtained.  In  particular, 
we  show  that  to  the  accuracy  of  the  calculation,  a  simple  (practically  constant) 
trial  function,  plus  one  iteration,  yields  a  decay  time  identical  with  the  re- 
s\ilt  obtained  by  Wrubel. 

We  note  that  this  report,  and  the  preceding  one,  may  be  read  indepen- 
dently. 
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2.  Assumptions 

In  order  to  simplify  the  calcvilation,  ve  shall  make  certain  assump- 
tions concerning  a)  magnitude,  and  b)  mathematical  form.  In  this  section  we 
discuss  our  assumptions,  their  justification,  and  the  sirr.plifications  that 
result. 

If  c  is  the  velocity  of  light,  X   the  decay  time,  C  the  conductivity 
(e.s.u. ),  and  L  a  typical  dimension  of  the  conducting  region, 

(2.1)  r  ^  ^- . 

c 

We  now  assume  that  the  decay  times  are  always  very  largej  that  is,  using  (2.1), 

■«- 
we  require 

(2.2)  (^  ^c-r»i  . 

The  first  form  states  that  during  the  decay  time  an  electromagnetic  disturbance 
can  travel  a  distance  much  larger  than  Lj  the  second  form  states  two  equivalent 
facts:  a)  the  conduction  current  is  much  larger  than  the  displacement  current, 
and  b)  the  volume  charge  decays  in  a  time  very  much  less  than  the  decay  time. 

In  Gaussian  units,  with  e  =  jj.  =  1,  the  magnetic  field  is  always 
divergencelessj  because  the  volume  charge  disappears  so  qiiickly,  the  electric 
field  may  also  be  regarded  as  divergenceless.  Thus, 
(2.3)  V  .  H  =  0  , 

{2,k)  V  .  E  =  0  . 

Since  the  fields  decay  as  e~    ,   the  first  Maxwell  equation  becomes: 

H 
(2.5)  ''<£  =  —   , 


In  oitr  problem,  using  the  sun's  radius  for  L,    \~ry     -^  10"^  •     ^^  *he  center 
of  the  sun,     (7''^'^  ICr   j   this  falls  to  ICr     on  the  shell  having  a  radius  90 /o 
of  that  of  the  sun. 
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In  the  second  Maxwell  equation,  we  omit  the  displacement  c\irrent: 

(2.6)  VX  H     =       lil   Q-E     , 

—  c  — 

E 
The  neglect  of  the  displacement  current,  — — •  ,  in  (2,6)  is  justified  by  (2,2): 

C  t 

according  to  (2,5),  E/n/  H  —p-     j  therefore,  in  (2,6),  the  displacement  current 
is  of  order  l"- rr}   compared  with  the  curl  H  terra,  and,  if  the  restriction 
(2,2)  applies,  may  be  dropped.  The  theory  resulting  from  Eqs,  (2.5)  and  (2,6) 
is  sometimes  called  the  'quasi-steady'  theory. 

Eqs,  (2,3)  and  (2,5)  are  consistent,  Comi^aring  Eqs,  (2,Li)  and  (2,6) 
leads  to  an  additional  relation: 

(2.7)  E  •  V  C  =  0  . 

In  oiir  problem,  O^  is  spherically  symmetric j  (2,7)  therefore  requires  that  E 
have  no  radial  conponent  (H-modes  only). 

Outside  the  conducting  region  (radius  R,  surface  S,  volume  V)  the 
above  equations  (with  C  =  O)  still  apply,  at  least  up  to  distance L  that  satisfy 
(2,2) J  beyond  these  distances,  the  field  is  negligibly  small.  That  is,  outside 
the  conducting  sphere  both  E  and  H  satisfy  Laplace's  equation,  with  the  two 
fields  coupled  according  to  (2,5) o 

The  E  field  in  the  external  region  may  be  replaced,  in  its  influence 
on  the  E  field  in  the  conducting  region,  by  a  boundary  condition.  We  express 
the  boiindary  condition  in  terms  of  the  unit  normal  n  (pointing  out  of  the  con- 
ductor) and  the  (symmetric)  dyadic  y  : 

(2.8)  (V  X  E)  X  n  =  y'.  E  (on  S), 

[31 
The  notation  has  been  chosen  to  agree  wiuth  that  used  previously"-  -' , 
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Eliminating  H  from  (2.5)  and  (2,6), 

(2.9)  (v^  ^  _Un^  E  ,  0^ 

where  the  Laplacian  operates  on  the  rectangular  components  of  E,  Eqs.  (2«)4), 
(2,7),  (2,8)  and  (2,!?)  constitute  a  vector  eigenvalue  problem.  We  now  reduce 
it  to  a  scalar  problem  by  specializing  the  form  of  E. 

Introdticing  a  spherical  coordinate  system  (r,©,0)  \dth  unit  basis 
vectors  e  ,  e^,  e^^  >  we  limit  the  electric  fields  to  the  forms  considered 

~r2i 

by  Wrubel  ^-^  ,   namely, 

(2.10)  E  =  E^   ^  (l^)  %f  • 

Here,  [i   =  cos  6,  and  P  is  the  associated  Legendre  ftmction  of  degree  n  and 
order  1. 

Fields  of  the  form  (2,10),  since  they  are  independent   of  0,   satisfy 
i2»h)i   since  they  teve  no  radial  component,  they  also  satisfy  (2,7),  All  that 
remains  is  to  satisfy  Eqs,  (2,8)  and  (2,9), 

Consider  (2,9).  In  terms  of  the  rectangular  xmit  vectors  e  ,  e  ,  e  , 

(2.10)  becomes: 

(2.11)  E  =  E^  Re|:cJ;  (©,!^(ey+  ^^Jl, 

I^ 

order  1, 

Using  (2,11)  in  (2,9): 
,2 


where  Re  means  'real  part  of  and  x  is  the  spherical  harmonic  of  degree  n  and 


(2,12)       14^  -  ^ii2^  p(r)  +  42- C^  (r)p(r)  =  0. 


dr        r"         c~r 
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Now  we  t\im  to  (2,8).     Outside  the  conductor,     C"=  0  in  (2.12), 
Therefore, 

(2,13)  p(r)  ^  \     i  r  >  R  . 

r 

Putting  (2.10)  in  (2.8),  the  boundary  caidition  reduces  to: 

y^^  is  evaluated  from  (2.13);  (2,lU)  then  becomes: 

(2.15)  I  §  +  n  -  0    J  r  .  R  . 

The  electric  field  must  be  finite  at  the  origin;  therefore, 

(2.16)  p(0)  =  0  . 

We  now  have  a  scalar  eigenvalue  problem:  namely,  to  find  those  values  of  f 
for  which  there  exists  a  solution  of  (2.12)  subject  to  the  boundary  conditions 
(2,15)  and  (2,16). 

r2i 

Wrubel'-  ■'  solved  this  problem  by  assuming  a  value  for  'T  in  (2.12) 
and  integrating  (2.12)  numerically  from  r  =  Otor=Rjp  and  t^  at  r  »  R 
were  then  tested  by  compaiT.son  with  (2,l5)j  if  (2,l5)  was  not  satisfied,  the 
process  was  repeated  with  a  new  value  assumed  for  f  until,  eventually,  all 
conditions  of  the  problem  were  satisfied, 

3«  Variational-iterational  procedure^ 

In  this  section  we  shall  set  up  a  variational  iterational  procedure 
for  solving  the  scalar  eigenvalue  problem  presented  by  Eqs.  (2.12),  (2.l5), 
and  (2,16)  for  n  =  1, 

'  The  general  method  is  given  by  G,  Ter\ple  and  W.G,  SLckleyL'^-l.  The  applica- 
tion to  magnetic  decay  problems  is  discussed  in  Ref.  [^].  Several  variational 
principles  are  given  in  Ref,  [3]}  in  particular,  a  (vector)  variational  itera- 
tional procedure  for  computing  t    from  the  electric  field  distribution  alone 
is  described;  if  (2J.0)  is  substituted  in  this  variational -iterational  scheme, 
the  scalar  equations  of  this  section  follow. 
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JFlrst  consider  the  variational  principle: 


(3.1)       ^ 


R    C  „  2       ,2  1 


c  f  >..        p 

dr  G-  p 


In  (3.1),  6f  =  0  for  variations  of  p  about  the  solution  of  (2.12)  subject  to 
the  restrictions  (2.15)  and  (2,16)^-''   (putting  n  -  l). 

The  ^ectrura  of  t   is  discrete  and  bounded  from  above  ^  -' .  We  limit 
ourselves  to  the  largest  eigenvalue  for  n  ■=  1  and  call  this  eigenvalue  T^  ^ 
(no  node).   The  value  of  f  obtained  by  inserting  a  trial  field    in  (3.1)  will 

^^  ^  '^1,0- 

Call  this  initial  trial  field  p- .  A  sequence  of  iterated  functions 

p,  (k  =  2,3,.»«)  is  defined  by  those  solutions  of 

(3.2) 

that  satisfy 
(3.3) 


(3.U) 


-\^.- 

c 

P^Co)  - 

0 

-1       5 

r  =  R  . 


'k   dr 


An  explicit  solution  of  these  equations  is  provided  by  the  Green's  function 


(3.5) 


2  , 

G(r,r  )   =  — f   J  r  <  r 

3r 

«2 

3r    *  r  >  r  . 


*•«■ 


Assumed  not  orthogonal  to  the  eigenfunction  belonging  to  ^^   q. 
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Since  G(r,r  )  satisfies  (3.3),  (3.)4)  and 

2 

(3.6)  -^  G(r,r')  -  -^  G(r,r')  =  -5(r-r') 
dr  r 

we  find 

R 

(3.7)  Pi^(r)  =   r  dr'  G(r,r')  ^  Cr  (r')  pj^_^(r')  . 

Jo  c 

According  to  (3.5),  G(r,r  )   =  G(r  ,r),     Tlierefore,  from  (3.7), 

R  R 

(3.8)  r    dr  pj^(r)  (7(r)  P^_3_(r)     =        f    dr'  P^_-^(^' )  C'(r')  p^(r')     =     (k+i^-l), 

where  the  last  notation  emphasizes  that  only  the  svm   of  the  subscripts  enters 
into  the  value  of  the  integral. 

We  put  p  =  p,  in  (3.1)  and  label  the  resulting  value  of  the  decay 
time  '^oi^.  Using  (3.2),   r 


[3]+ 


It  may  now  be  shown '--J  ^  that,  as  k  ->oo,  the  11,  increase  in  value, 
approaching  T  ^  as  a  limit,  while  the  functions  p,  converge  towards  the  exact 
eigenfunction  corresponding  to  T  „. 
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li.     Numerical  restilts 

We  have  already  taken  n  =  1,  If  we  introduce  a  dimensional  variable 
X  by 

P. 

the  integrals  in  (3»6)  become: 


(U.l)  x=|;  0<x<l, 


a.2)        (k*^)  =  n7{0)     j  dx  -J[^  Pj^(x)5^(x)  . 

The  variation  of  the  function   ■  >^y  in  the  sun  (  in  Wrubel's  notation  0  '  ) 
has  been  tabiilated  by  Wrubel  in  steps  of  .1.  The  values  given  by  Wrubel  are 
reproduced  in  'Eable  I.  Since   -/^^  is  only  available  numerically,  the  in- 
tegrals of  the  form  of  (ii.2)  that  occur  have  been  evaluated  by  numerical  methods. 
The  procedure  has  been  to  apply  Simpson's  Rule,  using  as  a  net  the  10  subdivisions 
appearing  in  Table  I.  Note  that  the  values  taken  on  by  the  p.  at  x  =  1  are  un- 
important  since  the  p,  occur  in  (U.2)  multiplied  by   g./^-,  and  this  factor 
vanishes  at  x  «  1.  The  practical  result  of  this  is  that  the  boxindary  condition 
(3»U)  may  be  ignored  in  this  problem. 

At  X  =  0,  however,  the  boundary  condition  must  be  included.  We  have 
done  this  for  p^ (x)  by  choosing  this  function  to  be: 

(h.3)  P^M 

1     i       X  /O 

Using  (U.2)  and  (li.3),  we  find: 
(il.U)  (2)  «=  .202  R  O-(O)   . 
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Iterating  p, (x)  by  use  of  (3»7),  vre  find: 


f  \    1  n  r2  'y(o)   ]l 


c 

(l4.5) 


i/„-^'^i{^v^''-^X-"^7^^^^'' 


Un  R    — ^ —     a^  W 
c 


thereby  defining  a  (x)  ♦  Values  of  a^Cx)  computed  from  (1^.3)  and  (h«5)  are 
listed  in  Table  I. 

With  Pp(x)  available,  we  obtain  frora  (Ii.2), 

([,.6)  (3)     =    1.087X10"^     l;n  R^    -2^         R  O-(O) 

c 

(I4.7)  (i4)     =    6.514   >i  10"^    Un  R^  -^1^  j  R   G-iO)   . 

Before  combining  these  values  according  to  (3»?)  we  must  have  a  value 

2       2  m 

for  the  factor  i^nR  O'(0)/c  .  In  Wrubel's  work  this  factor  (PTin  Wrubel's 

notation)  has  the  value   (see  Wrubel's  Table  3) 

(ii.S)  Un  R^    -^1^      =       6.56  X  lO"'-'-  years 

c 

2 

(Note  that  Wrubel's   CT"  and  ours  differ  by  a  factor  c  ). 

Collecting  these  results,  we  obtain  two  estiinates  for    T,   ^z 
(li.9)  T    =    3.53  X  10^    years 

(I4.IC)  tj^  =    3.96  X  10^    years     . 
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As  expected,  T,  >  X«  Wrubel's  value  for  the  decay  time  is: 
(li.ll)  '^Wrubel  ^  ^'^^  ^  ^^^   y^^^^  • 

This  differs  by  3  /^  from  f,  . 

However,  we  should  note  that  Wrubel  uses  a  subdivision  that  employs 
10  times  the  nvimber  of  points  we  use.  To  evalioate  the  effect  of  our  coarser 
net,  we  have  used  Wrubel's  'exact'  values  for  p(x)  and  ^  in  (3«l)«  As  in 
the  previous  evaluations,  the  integrals  were  done  by  Simpson's  Rule  i:isLng  only 
the  ten  subdivisions  given  in  Table  I,  Since  it  is  this  value  of  f ,  rather 
than  (lj..]l),  that  should  be  compared  with  the  variational  results  t^y    ^j  , 
we  call  this  decay  time  T    .  •  We  find: 

(1).12)  z^y,^^^     =  3.96X10^  years  . 

Thus,  the  3%  difference  between  (h.ll)  and  (h.lO)  may  be  considered 
as  arising  solely  from  the  fact  that  the  net  used  for  (It. 11)  was  10  times  finer 
than  that  used  for  (U«10) . 

Vfe  conclude  that  the  use  of  the  simple  trial  function  (l4«3),  plus  one 
iteration,  in  the  variational  principle  duplicates,  to  the  accuracy  of  the  calcu- 
lation, the  value  Wrubel  obtained  by  numerical  integration  of  the  differential 
equation. 
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TABLE  I 


CM 


OTO)  "2 


aAx) 


0.0  1.000  0.00 

.1  .808  3.99  X  10"^ 

.2  .li7U  6.55  X  10"^ 

.3  .26U  7.60  X  lo"^ 

.li  .1U6  7.70  X 10"^ 

•5  .0816  7.3li  X  10"^ 

.6  .cUli7  6.7U  X  lo"-^ 

.7  .0231  6.07  X  10"^ 

.8  .0103  5.Uh  X  lo"^ 

.9  .00305  li.6l4  X  lo"-^ 

l.C  .00000  U,36  X  lo"^ 
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